Additional is-wave Cooper pairing driven by nematic order in d-wave superconductors 
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We make a non-perturbative analysis of the strong interaction between gapless nodal fermions and 
nematic order parameter in two-dimensional d x 2_ y 2 superconductors. We predict that the critical 
nematic fluctuation can generate a dynamical nodal gap if the fermion flavor N is smaller than a 
threshold N c . Such gap generation leads to an additional is-wave Cooper pairing instability, and 
induces a fully gapped d x 2_ y 2 +is superconducting dome in the vicinity of nematic quantum critical 
point. The opening of dynamical gap has important consequences, including saturation of fermion 
velocity renormalization, weak confinement of fermions, and suppression of observable quantities. 
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One of the most prominent properties of high-T c 
copper-oxide superconductors is that they exhibit a num- 
ber of long-range orders upon changing the chemical 
doping, such as antiferromagnetism, superconductivity, 
stripe, nematic state, and so on. The competition and 
coexistence between superconductivity and other long- 
range orders are believed to be fundamental issues in the 
studies of high-T c superconductors. Among the orders 
that are in competition with superconductivity, the ne- 
matic order, which spontaneously breaks C4 symmetry 
down to C 2 symmetry, has attracted special theoretical 
and experimental interest in the past decade 

In recent years, strong anisotropy in electronic proper- 
ties has been observed in various experiments performed 
on YBa 2 Cu 3 6 +5 [H3 and Bi 2 Sr 2 CaCu 2 8+(5 Q. Such 
anisotropy is widely attributed to the formation of an 
electronic nematic state. An important indication of 
these experiments is the existence of nematic quantum 
phase transition in these two typical superconductors. 
Interestingly, similar nematic states are also observed in 
several other correlated electron systems, including iron- 
based superconductor Q , heavy fermion superconductor 
P3, Sr 3 Ru 2 7 superconductor (Tlj . and semiconductor 
heterostructure [12 ], 

High-T c superconductors are known to have a d x 2_ y 2 
energy gap, which vanishes linearly at four nodal points 
(±^,±-|). Due to this property, gapless nodal quasi- 
particles (qps) are present in the superconducting dome 
even in the low-energy regime. When a nematic quantum 
phase transition occurs in the superconducting dome, the 
fluctuation of nematic order parameter couples to the 
gapless nodal qps. This coupling becomes singular at the 
nematic quantum critical point (QCP), and can lead to 
unusual behaviors. Vojta et at analyzed this model by 
means of e = 3 — d expansion and found runaway behav- 
13j. Later, perturbative expansion in terms of 1/N 



lor 

with N being the fermion flavor has been extensively ap- 
plied to address this issue [I3 - l20| . For instance, Kim 
et. al. revealed a second-order nematic phase transition 
after performing a large- TV analysis More recently, 
renormalization group calculations of Huh and Sachdev 
[l5l | discovered a novel fixed point with extreme fermion 
velocity anisotropy, i.e., S = v^/vp — > 0, where Vf,a are 



Fermi/gap velocities of nodal qps respectively. Subse- 
quent studies showed that such extreme anisotropy can 
give rise to a variety of nontrivial properties, such as 
non- Fermi liquid behaviors [l6| , enhancement of thermal 
conductivity 17 1, suppression of superconductivity [l9| . 
and first order instability [20j . 

We should note that all pervious field-theoretic analy- 
sis are based on perturbative expansions [l3Tj20| . In this 
formalism, one starts from the free fermion propagator, 
G ! Q 1 (w,k) = —ioj + VFk x T z -\-VAk y T x , and then calculates 
the one- loop fermion self-energy S(w,k) that is induced 
by the interaction with nematic fluctuation. Generically, 
S(u;,k) can be expanded as 

k) = — zSow + Yi\Vpk x T z + £ 2 ?j Ak y T x '', (1) 

where the functions So. 1,2 ar e the temporal and spatial 
components respectively. The fermion damping effect is 
encoded in So [3 , whereas the velocity renormalization 
can be obtained from £0,1,2 T5|. However, in principle 
there could be the fourth term, mr y , corresponding to 
a nonzero mass gap term of nodal qps. This mass term 
can never be obtained to any finite order of perturba- 
tive expansion, but may be dynamically generated if we 
perform non-perturbative calculations. 

In this paper, we study dynamical gap generation of 
originally gapless nodal qps due to nematic fluctuation 
by means of non-perturbative expansion. With the help 
of Dyson- Schwinger (DS) equation that connects the free 
and complete propagators of nodal qps, we obtain a non- 
linear gap equation of fermion mass m in the vicinity of 
nematic QCP. After solving this equation, we find that a 
nonzero mass gap, mr y , is dynamically generated when 
the fermion flavor N is below certain critical value N c , 
i.e. N < N c . We demonstrate that the dynamical gap 
to induced by nematic order corresponds to a secondary 
is-wave Cooper pairing formation, so the critical nematic 
fluctuation drives a transition from a pure d x 2_ y 2 super- 
conducting state to a d x 2_ y 2+is superconducting state in 
the vicinity of nematic critical point. As a consequence, 
the superconducting state is fully gapped and the massive 
nodal qps are weakly confined by a logarithmic poten- 
tial. Moreover, the dynamical gap leads to saturation of 
fermion velocity renormalization and strong suppression 
of some observable quantities. 
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We start from the following action [14H1 

S = + S$ + s. 

where the free action for nodal qps is 

-*1 q (-kj + v F k x T z + v A k y T x )^ la 



(2) 



S, 



(2tt) 2 2tt 
d 2 k duj 



(2tt) 2 2tt 



$i a (-W + V F k y T Z + V A k x T X )^2a(3) 



Here, the Nambu spinor \fr is defined as, "Fia = 
(/ia, £abfl b ) T and * 2a = (/2a, fab/],,) 7 , where e ab = -e ba 
with spin indices a, 6. The operators /1, /2, /3 and /4 are 
nodal qps excited from four nodal points (f , f ), (—§,§), 
(— f ) — f ) ; and (§,—§), respectively. The free action for 
nematic order parameter cf> is 



d 2 xdr 



4! 



where </> is a real scalar held since the nematic transition 
is accompanied by a discrete symmetry breaking (i.e., 
Ising-type) . The coupling between the nodal qps and the 
nematic order parameter is [l3| 

= J d 2 KdT{\^\ a T X ^ la + *L^*2a)}- (5) 

The physical flavor of nodal qps is determined by 
spin degeneracy, so W = 2. However, in order 
to carry out analytical calculations, it is convenient 
to assume a large TV and then perform 1/N expan- 
sion. The free propagators for nodal qps and 
are G% (w, k) = {—iu + VFk x T Zj rv A k y T x )~ 1 and 

G% (w,k) = (—iuj + vpkyT 7 - + v A k x r x ) 1 , respectively. 
To the leading order of 1/N expansion, the polarization 
function due to particle-hole excitations is given by 



11(0, q) = A 2 EE /(0^^>< k ) 



a=li=l,2" v ' 

After straightforward calculations, it is easy to get 
NX 2 n 2 +v%ql 



(6) 



n(o, q ) 



16vfva 



<li 



+ {qx^q y ), (7) 



where qi = y ft 2 + v 2 F q\ + v\q 2 . The polarization is lin- 
ear in \q\, and dominates over the kinetic term q 2 at low 
energies. We can drop the q 2 term and write the effective 
nematic propagator as 



rr 1 (n,q) = r + n(n, q ). 



(8) 



The free fermion propagator will receive corrections 
due to fluctuation of nematic order parameter. After 



including these corrections, the complete propagator of 
takes the following general form 



G^/ la (u>, k) 



1 



-iuiA + vpk x AiT z + v l \k v A2T x + mTV ' 



(9) 



where ^0,1,2 are wave- function renormalizations and m = 
m{uj,k x ,ky) is a fermion gap. The mass gap term can 
never be generated as long as the fermion self-energy is 
calculated perturbatively by using the free fermion prop- 
agator G°(w, k). To examine the possibility of dynamical 
gap generation, we should instead utilize the following 
self-consistent DS equation 

G-l(e,p)= [G^ la (e,p)] _1 -E la (e,p) (10) 

where the self-energy is computed using the complete 
fermion propagator, 



£ia(e,p) 



1 



n(e u), p k) 



(11) 



Substituting the complete propagator Eq.(9) into the DS 
equation, one can derive four self-consistently coupled 
equations of Ao,i,2 and m. Compared with the equation 
of m, the equations of ^0,1,2 contain an extra pre-factor 
of 1/N, and are therefore relatively unimportant at large 
N. To the leading order of 1/N expansion, we assume 
that Aq — A\ — A2 = 1 , and then obtain a gap equation, 



= A 



m(s,p x ,p v ) 
duj d 2 k 



m(uj,k x ,ky 



2tt (2tt) 2 cj 2 + v 2 F k 2 x + v 2 A kl + m 2 {uj, k x , k y ) 
1 



r + n(e — u>, p k) 



(12) 



If this equation has only vanishing solution, m = 0, then 
the nematic fluctuation can not open any gap. A fermion 
gap is dynamically generated once this equation develops 
a nontrivial solution, m ^ 0. Due to the anisotropic 
nature of nematic fluctuation, the integrations over uj, k x , 
and k v have to be performed separately, which greatly 
increases the difficulty of numerical computations. For 
simplicity, here we consider the isotropic limit, Vf = v&. 
In this case, the dynamical gap becomes m(e, |p|), and 
the polarization is simplified to 



n(n,q) 



AA 2 2ft 2 



4|q| 2 



16v 2 F 



(13) 



We first consider the nematic QCP and take r = 0, 
then the pre-factor A 2 on the right-hand side of Eq.(12) 
cancels exactly the factor A 2 appearing in polarization 
n(e — u>, p — k). As a result, the gap equation is actu- 
ally independent on coupling constant A, and the critical 
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FIG. 1: (a) Relationship between dynamical gap m(0, 0) and 
fermion flavor N at nematic QCP, r = 0; (b) Dependence of 
m(0, 0) on tuning parameter r for N = 2. Dynamical gap is 
destroyed once r exceeds certain critical value r c . 



point of dynamical gap generation is solely determined 
by flavor N. After numerical computations, we find a 
critical fermion flavor N c » 2.4 at r = 0. The nodal 
qps remain gapless, m = 0, when iV > iV c , but acquire 
a nonzero dynamical gap m when A < A c . Figdfa) 
presents the dependence of gap m(v = 0, k = 0) on fla- 
vor N. It is clear that the dynamical gap decreases very 
rapidly as flavor N grows, and vanishes continuously as 
N — > N c . It is also interesting to examine how the tuning 
parameter r affects the gap generation. Actually, if we 
move away from the nematic QCP, r becomes finite and 
the nematic fluctuation is no longer critical. For r ^ 0, 
the coupling parameter A cannot be exactly canceled, but 
one can absorb it into r by taking r — > r/A 2 . We find that 
the dynamical gap is significantly suppressed by growing 
r, and completely destroyed once r exceeds certain crit- 
ical value r c , which is shown in FigQJb). We therefore 
conclude that the dynamical gap generation is mediated 
by the critical fluctuation of nematic order parameter, 
and exists only in the vicinity of the nematic QCP. 

Now let us discuss the physical meaning of the dynam- 
ical gap. The nonzero gap m of nodal qps generates an 
extra term which should be added to the Hamiltonian, 



Hm = 



d 2 k 

(27T) 2 

W) 



2 { lTn [(/aVit + foht 



(fyit + fuf*t)+0-** 2,3^4)}. 



One can see that such term corresponds to the formation 
of singlet Cooper pairs between the gapless qps excited 
from opposite nodes. We therefore have obtained a sec- 
ondary, nematic-order driven, is-wave superconducting 
instability on top of the pure d x 2_ y 2 superconductivity. 
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FIG. 2: Schematic phase diagram. x c denotes the nematic 
QCP. The shadowed region around x c represents the emer- 
gent fully-gapped d x 2_ y 2 +is superconducting dome within a 
larger, pure d x 2_ y 2 superconducting dome. 



As already pointed out, this dynamical gap is opened 
only when r is zero or very small, and is rapidly destroyed 
as r increases, which implies that the secondary is- wave 
superconductivity is achieved only in the close vicinity of 
nematic QCP. Upon approaching the nematic QCP, there 
is a zero-temperature phase transition from a pure d x 2_ y 2 
superconducting state to a fully gapped d x 2_ y 2 + is su- 
perconducting state. At finite temperature, T / 0, the 
critical nematic fluctuation is weakened strongly due to 
the thermal screening effects, hence the dynamical nodal 
gap cannot survive at high temperatures. According to 
these analysis, we can infer that a small d x 2_ y 2 + is su- 
perconducting dome emerges around the nematic QCP, 
which is schematically shown in Fig. 2. It is interesting 
to notice that such nematic fluctuation-driven supercon- 
ducting dome is analogous to that is formed on the border 
of an antiferromagnetic quantum critical point in some 
heavy fermion compounds [2lj . 

Now we discuss the consequences of nonzero dynamical 
gap. For simplicity, let us assume a constant gap m, 
which gives rise to a new polarization function 



n(q,fi) 



nx 2 n 2 



2ttvfva 




1 m 




m 2 \ 




G- 





In the low 
n(q,Q) 



energy 

NX 2 



limit, 
■(2ft 2 - 



it 



+ (<7x ** %)■ 



takes the simplified form 
Using this polariza- 



v 2 M 2 



tion, it is easy to obtain the following effective potential 
[22| between two massive nodal qps, 



V(R) oc 



d 2 q e 



q-R 



6ua 



(2tt) 2 n(q) NX 2 v F 



mln(mR). (14) 



This potential grows logarithmically as the distance 
R between two fermions is increasing, so the gapped 
fermions are weakly confined. Such a gap-induced con- 
finement of fermions is similar to that in QED3 [22| ■ 
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FIG. 3: Momentum dependence of fermion velocities. Blue 
solid line: 8=1 and m = 0; Red Solid line: 8 = 0.1 and 
m — 0; Blue dashed line: 8 = 1 and m/v^A = 10~ 6 ; Red 
dashed line: 8 = 0.1 and m/iifA = 10 -6 . 



When the nodal qps are massless, their velocities are 
renormalized by the nematic fluctuation and thus become 
strongly momentum-dependent, Vf,a — ► VF,A(k). Both 
of «f,a(&) vanish as k — > 0, but v&(k) vanishes much 
faster than vp(k) (23j |. leading to extreme anisotropy, i.e., 
S(k) = VA(k)/v F (k) -> [l5l - [T9| . Nevertheless, once a 
finite gap m is opened, there is no longer strong renor- 
malization of velocities. As shown in Fig. 3, the velocities 
vf.a are still fc-dependent and decrease with decreasing 
k at high energies, but saturate to finite values once k is 
smaller than the scale set by to. Therefore, the singular 
velocity renormalization and the extreme anisotropy are 
both prevented by the dynamical nodal gap m. 

The dynamical gap also has important influence on 
observable quantities. For example, the density of state 



of nodal qps is linear in energy, p(oj) oc |w| when to = 0, 
but becomes p(cj) cx |a>|0(|o;| — to) when to ^ 0, which 
vanishes for |w| < to. Accordingly, the specific heat is 
strongly suppressed as Cy oc m 4 /T 2 exp(— m/T) in the 
temperature T <C to. In addition, for finite to. the low 
temperature thermal conductivity is of the form [2J] , ^ = 

%"(^f + l^) r'i+mi > wnere r is impurity scattering rate. 
If m = 0, thermal conductivity is universal and impurity 

independent, i.e., f = ^(Zz + m.) [2^|. In contrast, for 
finite to, thermal conductivity is suppressed. 

The nematic state is indeed equivalent to a supercon- 
ducting state with d x 2_ y 2 + s gap, which was pointed 
out in Ref. [l3| . Therefore, in the vicinity of a quantum 
critical point between a pure d x i_ y i + s superconducting 
state and a d x 2_ y 2 +s superconducting state, the singular 
fluctuation of s-wave order parameter can also leads to a 
fully gapped d x 2_ y 2 + is superconducting state, provided 
that fermion flavor N is sufficiently small. 

According to our results, it turns out that the fermion 
flavor N is a crucial parameter which determines the low- 
energy behaviors caused by the nematic order. When 
N > N c , the non-perturbative effects of nematic fluctua- 
tion is unimportant, so one can trust the results obtained 
by perturbative calculations, such as extreme anisotropy 
and other unusual properties [H ll6l - [2f]| . If TV < N c , 
however, the non-perturbative effects become significant, 
and can drive additional is- wave superconducting pairing 
between the originally gapless nodal qps. Our leading- 
order calculations found that N c ps 2.4, which is larger 
than physical flavor N — 2. It should be interesting 
to study how N c is quantitatively affected by high or- 
der corrections in the future. In any case, our prediction 
of a nematic order- induced d x 2_ y 2 + is superconducting 
dome is novel and would shed light on the investigation 
of nematic order in correlated electron systems. 
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